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METHOD OF ANALYZING CORROSION AND CORROSION PREVENTION 



Technical Field 
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The present invention relates to an analyzing 



method, which uses a computer for predicting corrosion and 
corrosion prevention. More particularly , the present 
invention provides an analyzing method that is preferably 
applicable to the problem of macro-cell corrosion (cathodic 

10 protection) such as heterogeneity metal contacting corrosion 
such as galvanic corrosion and differential aeration 
corrosion among the problems of metal corrosion and 
corrosion prevention. The present invention is also 
applicable to a system in plating, battery cell, and 

15 electrolytic tank, where an anode and a cathode are 
macroscopically disposed with an electrolyte interposed 
therebetween to form a potential field, other than the 
problems of metal corrosion and corrosion prevention. 

20 Background Art 



such as seawater, an object tends to suffer damage due to 
macrocell corrosion such as galvanic corrosion caused when 
different metal materials are used or flow velocity 
25 corrosion caused by an irregular flow velocity distribution 
(differential aeration corrosion caused by a flow velocity 
difference) . It has been desired to accurately predict 
those types of corrosion in advance and provide 



In a solution having a high electric conductivity 
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countermeasures against them. ^Cathodic protection'' based 
on the positive utilization of a corrosion inhibiting 
phenomenon on the cathode of a macrocell is widely used as a 
most basic corrosion prevention process. It is required for 
5 cathodic protection to predict a corrosion prevention range 
and the rate of consumption of a sacrificial anode depending 
on the material and location of the anode, the shape, 
material , and arrangement of a device to be protected 
against corrosion, and solution conditions (electric 

10 conductivity and flow velocity) . 

There are limitations on an experimental approach 
to the prediction of a macrocell because the shape of the 
field has a large effect on the behavior of the macrocell. 
Specifically, for example, even if an experiment is 

15 conducted on galvanic corrosion and the effects of various 
factors including an area ratio, a combination of materials, 
and the electric conductivity of a solution are studied in 
detail, the results are only applicable to the three- 
dimensional shape of a region that is occupied by the 

20 solution in the experiment. Since actual devices and 
structures are of complex shapes, it is not possible to 
accurately estimate a liquid junction resistance in a 
macrocell, and it is difficult to apply experimental results 
directly. Furthermore, each time the shape of a device to 

25 be protected against corrosion is changed, it is practically 
impossible to conduct an experiment by assuming the changed 
shape . 

Therefore, it has been customary to predict 
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macrocell corrosion and cathodic protection on a real 
structure according to empirical rules in most cases. Many 
attempts have been made to perform more accurate, 
quantitative predictions. According to one effort, 

potential and current density distributions are determined 
by purely mathematically solving a Laplace equation, which 
governs a potential distribution. However, objects to be 
analyzed are limited to relatively simple systems such as 
flat plates, cylinders, etc. While a conformal mapping 
process and a process using an electrically conductive paper 
have been employed for a long time as a process for 
analyzing an electric field problem, these processes can 
handle only a two-dimensional field. 

With the development of the computer technology in 
15 recent years, attempts have widely been made to apply a 
numerical analysis based on a differential method, a finite 
element method, and a boundary element method. However, the 
differential method and the finite element method are 
disadvantageous in that they require a long computational 
20 time as the overall object has to be divided into elements. 
According to the boundary element method, on the other hand, 
the time required for dividing the object into elements and 
performing computations can greatly be shortened because 
only the surface of the object needs to be divided into 
25 elements. Considering that the boundary element method is 
an optimum process for analyzing a corrosion problem where 
physical quantities on the surface such as a potential and a 
current density are important, the inventor of the present 
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invention has developed an analyzing technique, which has 
applied the boundary element method to the prediction of 
macrocell corrosion and cathodic protection problems. 
[Fundamental equations and boundary conditions] 
5 The corrosion of a metal in an aqueous solution 

progresses as an electrochemical reaction in the form of an 
anodic reaction and a cathodic reaction, which are paired 
with each other. For example, with respect to the corrosion 
of iron in a neutral aqueous solution of salt containing 

10 dissolved oxygen, such as seawater, the reactions progress 
as indicated by the equations (1) and (2) : 

Fe -> Fe 2 + 2e~ (anodic reaction) (1) 
1/2 -0 2 + H 2 0 + 2e~ 20H" (cathodic reaction) (2) 
On a surface of metal, a region where an anodic 

15 reaction is taking place is referred to as an anode, and a 
region where a cathodic reaction is taking place is referred 
to as a cathode. Usually an anode and a cathode are minute 
and mixed with each other, and their positions are not 
constant. Therefore, corrosion progresses substantially 

20 uniformly as a whole while being subjected to some 
irregularities. However, if the material, the surface 
state, and the environment are not uniform, then the anode 
and the cathode are localized, and corrosion is concentrated 
on a certain location (anode region) . The former is 

25 referred to as microcell corrosion (cell refers to a battery 
cell) , and the latter is referred to as macrocell corrosion. 
Seawater pumps are often significantly damaged primarily by 
macrocell corrosion such as galvanic corrosion or 
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differential aeration corrosion. The cathode in macrocell 
corrosion is inhibited from corrosion as only a cathode 
current flows through the cathode. A method of corrosion 
prevention which positively utilizes this corrosion 
inhibiting phenomenon is cathodic protection. 

Either of systems of macrocell corrosion and 
cathodic protection can be considered as a cell in which an 
anode and a cathode are present with an electrolyte 
interposed therebetween. A potential «|>) in the electrolyte 
is governed by the following Laplace equation (3) : 

V 2 4> = 0 (3) 

It is assumed that, as shown in FIG. 1 of the 
accompanying drawings, the electrolyte is surrounded by 
boundaries T lr T 2r T 3a , and T 3c where T x represents a boundary 
(a boundary with a constant potential) where the value of 
potential <|> is fixed to (|> 0 , r 2 represents a boundary (a 
boundary with a constant current density) where the value of 
current density q is fixed to q 0 , and r 3a and T 3c represent 
the surfaces of the anode and the cathode, respectively. 

Boundary conditions at the respective boundaries 
are given by the following equations: 

On r\: O = 4> 0 (4) 

On T 2 : q {= Kd4>/3n} = q 0 (5) 

On T 3a : 4> = - f a (q) (6) 

On r 3c : 0> = - f c (q) (7) 
where K represents the electric conductivity of the 
electrolyte, d/dn represents a differential in the direction 
of an outward normal, and f a (q) and f c (q) represent 
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nonlinear functions indicative of polarization 

characteristics of the anode and the cathode, which can be 
determined experimentally. When the equation (3) is solved 
with respect to the equations ( 4 ) through ( 7 ) representing 
boundary conditions, potential and current density 
distributions in the vicinity of the surface can be 
determined. The potential § and an electrode potential E 
that is to be actually measured are related to each other by 
<J> = - E. 

[Solution according to the boundary element method] 

Upon ordinary formulation of the boundary element 
method, a boundary integral equation is derived according to 
the equation (3) : 

ck<j> = jfqdT- j<t>q*dr 

r r 

(8) 

where <|>* represents the fundamental solution of a three- 
dimensional Laplace equation, q* = K3<|>*/3n, T represents 
boundaries (= I\ + T 2 + T 3a + r 3c ) surrounding the 
electrolyte, and c is c = 1/2 for smooth boundaries and c = 
C0/27C at an angular point of angle CO. 

In order to numerically solve the boundary integral 
equation, it is necessary to make it discrete. When the 
boundaries are divided into many elements and ([>, q are 
approximated by discrete values at respective nodal points 
and interpolating functions, the following simultaneous 
system of algebraic equations is derived: 



bj 

/,<«,) 



(9) 



where bj (j — 1, 2, • • ■ , p) represents the value of a known 
component of or q on the boundaries r\ + Y 2r Xj (j = 1, 2, 
• * • , p) represents an unknown quantity corresponding to bj , 
5 fj (qj) (j =1, 2, s) represents a nonlinear function 

indicative of polarization characteristics , p and s 
represent the numbers of elements on the boundaries Ti + T 2 
and r 3a + r3 C , and [A] and [B] represent matrixes determined 
by the geometrical shape of the boundary T. Since this 
10 equation is nonlinear, repeated calculations are required to 
solve the equation. The inventor employs the Newton-Raphson 
method. 

[Process of analyzing axially symmetric region] 



15 components include axially symmetric regions, and it is 
desired to analyze these axially symmetric regions more 
simply. The following two methods are chiefly considered as 
a process of solving an axially symmetric problem. Namely, 
(i) a method of using a fundamental solution for an axially 

20 symmetric problem, and (ii) a method of using an ordinary 
fundamental solution for a three-dimensional problem and 
reducing the number of elements in view of axial symmetry 
when making boundaries discrete. Using a fundamental 
solution which meets an axially symmetric condition poses a 

25 problem in that integral calculations become more complex 
than using an ordinary fundamental solution. The present 
program employs a method of reducing the number of elements 



Many actual devices such as some pipe or pump 
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in view of axial symmetry when making boundaries discrete. 
This method will be described below. 



necessary to divide all boundaries into elements in order to 
make the boundary integral eguation (8) discrete. However, 
since <j> and q have the same value in the circumferential 
direction due to axial symmetry, the equation (8) can be 
modified as follows: 



where Ti D represents a range on a one-dimensional line. 
From the equation (10) , a simultaneous system of algebraic 
equations can be obtained simply by making Tid discrete. 
Therefore, based on axial symmetry, the number of unknowns 
can greatly be reduced, and the accuracy can be increased. 
[Region dividing method] 

For the sake of brevity, a region made up of two 
parts is considered as shown in FIG. 2 of the accompanying 
drawings. If an internal boundary plane is represented by 
r B , then since the equation (9) is satisfied in each of the 
regions, the following equation is satisfied: 
Region I 



In an ordinary three-dimensional analysis, it is 




(10) 




(11) 



Region II 
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(12) 

where the suffixes I, II represent quantities relative to 
the respective regions I, II, and the suffix B represents a 
quantity relative to the internal boundary plane T B . {X M } 
(M = I, II) represents a vector having as its elements 
quantities relative to boundaries other than the internal 
boundary plane T B among x ± and q if {b M } (M = I, II) 
represents a vector having as its elements known quantities 
corresponding to X M (or a function representing a 
polarization curve) . 

Since a potential and a current density are 
continuous in the internal boundary, the following equations 
are satisfied: 

4> IB = <f» IIB (13) 
q IB = - q IIB (14) 

In the equations (11) and (12), if [H m ] {cjP} (M = I, 
II) on the right side is transferred to the left side, and 
the equations (13) and (14) are substituted, then the 
following equations are obtained: 



[A'G IB -H 1B t 



f*1 



(15) 



(16) 



These equations can be combined into the following 
equation : 



A'G 18 -H IB Q 
0-G IIB -H m A n 





f > 






X 1 












~B l 0~ 








OB 11 




r 






X n 





b n 



(17) 

This equation is a nonlinear equation as with the 
equation (9) . 

As described above, the inventor has developed six 
programs for analyzing an open region (a region surrounded 
by an electrolyte that extends infinitely far away, as the 
outer surface of a ship) and a closed region (a region 
surrounding an electrolyte, as the inner surface of a pump) 
with respect to each of two-dimensional, three-dimensional, 
and axially symmetric objects, and solved practical 
corrosion and corrosion prevention problems. 

In actual system, some of six regions that can be 
modeled as two-dimensional objects (open and closed 
regions) , three-dimensional objects (open and closed 
regions) , and axially symmetric objects (open and closed 
regions) are present continuously. FIG . 3 of the 

accompanying drawings show a specific example. FIG. 3 shows 
a seawater pump of stainless steel having Zn sacrificial 
anodes disposed in three locations on the inner surface 
thereof and four Zn sacrificial anodes in the form of prisms 
spaced at equal intervals on the outer surface thereof. The 
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inner and outer surfaces of the pump are connected to each 
other via seawater, so that the inner surface should 
electrochemically affect the outer surface and the outer 
surface should electrochemically affect the inner surface. 
5 However, the seawater surrounding the outer surface 

of the pump occupies a wide region, and boundaries for 
dividing the outer surface are too large to handle it as a 
closed region. It is practically impossible to model and 
simultaneously analyze the outer surface as a three- 

10 dimensional closed region as with the inner surface. 
Therefore, the inner surface of the pump is analyzed as a 
three-dimensional closed region, and the outer surface of 
the pump as an open region. The inner surface of a guide 
casing is divided into seven flow paths by seven helical 

15 guide vanes. Since these flow paths are symmetrical, one of 
them is taken and divided into three-dimensional elements. 
The prismatic anodes on the outer surface of the pump are 
handled as axially symmetric on the assumption that strip- 
shaped anodes of the same area are mounted on the outer 

20 surface of the pump, and handled as open-region axially 
symmetric models. 

Actually, since the inner and outer surfaces of the 
pump electrochemically affect each other, they need to be 
analyzed taking that consideration into account. However, 

25 because analyzing programs for handling the respective 
regions are different from each other (the inner surface of 
the pump: a three-dimensional closed-region program, the 
outer surface of the pump; an axially symmetric open-region 
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program) , it has heretofore been impossible to perform an 
analysis taking the effect of the regions on each other into 
account. While the inventor has developed the region 
dividing method, as described above, the developed method 
5 can only analyze regions modeled in the same manner. 

Disclosure of Invention 

The present invention has been made in view of the 
above drawbacks. It is an object of the present invention 

10 to provide a method of analyzing corrosion and corrosion 
prevention in a system where two or more same or different 
regions of six regions that are modeled as two-dimensional 
objects (open and closed regions) , three-dimensional objects 
(open and closed regions) , and axially symmetric objects 

15 (open and closed regions) are present continuously, the 
method being capable of analyzing the two or more different 
regions that are present continuously. 

According to an invention defined in claim 1 , there 
is provided a method of analyzing corrosion and corrosion 

2 0 prevention in a system where two or more same or different 
regions of six regions that are modeled as two-dimensional 
objects (open and closed regions), three-dimensional objects 
(open and closed regions) , and axially symmetric objects 
(open and closed regions) are present continuously, 

25 characterized by the steps of dividing the system into 
regions and dividing each of the regions into elements 
depending on respective models (two-dimensional, three- 
dimensional, and axially symmetric) thereof, and using one 
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of the regions as an attentional region and another as a 
non-attentional region , and giving varying values of a 
current density or a potential to the elements where current 
densities and potentials on the elements which are 
5 positionally the same as each other on a division plane 
disposed between and shared by two regions are equivalent to 
each other , the relationship between the current densities 
and potentials on the elements being unknown, solving a 
discrete boundary integral equation corresponding to the 
10 non-attentional region for the given values of the current 
density or the potential thereby to express the unknown 
relationship between the current densities and potentials on 
the division plane with the known relationship between 
current densities and potentials on elements in the non- 
15 attentional region other than the division plane , performing 
a boundary element analysis on the attentional region using 
the determined relationship between the current densities 
and potentials on the division plane as a boundary condition 
on the division plane thereby to determine potential and 
20 current density distributions in the attentional region in 
its entirety f and performing a boundary element analysis on 
the non-attentional region again using the determined 
potentials or current densities on the division plane 
thereby to analyze the regions in a related manner. 
25 According to an invention defined in claim 2, there 

is provided a method of analyzing corrosion and corrosion 
prevention in a system where two or more non-attentional 
regions are present continuously with respect to the one 
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attentional region, characterized by the steps of 
determining the relationship between current densities and 
potentials with respect to both division planes in the 
method according to claim 1, analyzing the attentional 
region using the determined relationship as a boundary 
condition to determine potential and current density 
distributions in the attentional region in its entirety, and 
performing a boundary element analysis on the non- 
attentional region again using the determined potentials or 
current densities on the division planes thereby to analyze 
the regions in a related manner. 

According to the present invention, the discrete 
boundary integral equation corresponding to the non- 
attentional region is solved to express the unknown 
15 relationship between the current densities and potentials on 
the division plane with the known relationship between 
current densities and potentials on elements in the non- 
attentional region other than the division plane. In this 
manner, the relationship between the current densities and 
20 potentials on the division plane in the non-attentional 
region can be obtained, and the attentional region can be 
analyzed from the obtained relationship on the division 
plane taking the non-attentional region into account. 
Therefore , the different regions can be continuously 
25 analyzed in a relatively short time. 

The method according to the present invention is 
also applicable to a system where two or more non- 
attentional regions are present continuously with respect to 
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one attentional region. 



Brief Description of Drawings 

FIG. 1 is a diagram illustrative of boundary 
conditions for determining a potential or current density 
distribution ; 

FIG. 2 is a diagram illustrative of the division of 

regions ; 

FIG. 3 is a view showing a structure of a seawater 
pump as an example of an object to be analyzed; 

FIG. 4 is a diagram showing two different regions 
that are present continuously ; 

FIGS. 5A and 5B are views showing the division of 
regions of the seawater pump; and 

FIGS. 6A through 6C are views showing the results 
of an analysis of the seawater pump, FIG. 6A showing the 
shape of the pump (segment positions) , FIG. 6B showing a 
potential distribution on the inner surface of the pump, and 
FIG. 6C showing a potential distribution on the outer 
surface of the pump. 

Best Mode for Carrying Out the Invention 

FIG. 4 shows a region Q that is divided into two 
different regions Q lr Q 2 that are present continuously. 
Boundary elements T lr T 2 of those regions Q lf Q 2 are partly 
known. T B represents a boundary element on the dividing 
plane, which has a potential and a current density that are 
common as viewed from the region Q 1 or the region Q 2 . 
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When an ordinary boundary element method is applied 
to an attentional region, the following discrete equation is 
obtained: 



(18) 

where <|> 2 and q 2 represent a potential and a current density 
on the boundary T 2 , ife 2 and qe 2 represent a potential and a 
current density on the boundary T B as seen from the region 
Q 2, and [H 2 ] and [G 2 ] represent coefficient matrixes 
obtained by applying the ordinary boundary element method. 
If the entire boundary surrounding the region Q 2 is made up 
of n elements, then it is assumed that the n elements 
include m elements on the boundary T 2 and one element on the 
boundary r B . 

When the unknown boundary nodal point quantity is 
15 transferred to the left side and the known boundary nodal 
point quantity is transferred to the right side in the 
equation (18) to change the equation (18) to a form Ax = b 
and the boundary condition on the boundary T 2 is 
substituted, the following equation is obtained: 
20 x represents a vector composed of an unknown 

potential or current density, b represents a constant-term 
vector after the known boundary condition is substituted, 
and A represents a coefficient matrix determined by H, G 
matrixes and the gradient of the polarization curve. 
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(19) 



where <|> 2 ? , q 2 ? represent an unknown potential and current 
density on the boundary T 2 , x represents an m + 21 row 
vector, vector b represents an n row vector, and matrix A 
represents a matrix having a size of nx (m + 21) . 
The matrix A is divided as follows: 



A = 



a 



ln+l 



a 
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(20) 



When <j) 2 ? and q 2 ? are erased from the equation (20) , 
the relationship between <|> 2 B and q 2 B according to the 
following equation is obtained: 



(21) 



The equation (21) is an equation representing the 
relationship between the boundary nodal point quantities (j) 2 B 
and q 2 B on the boundary T B . This equation represents the 
relationship between (|> 2 B and q 2 B taking into account the 
effect of the non-attentional region Q 2 , and is considered 
as a boundary condition equivalent to the effect of the non- 
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attentional region. If this boundary condition is given as 
a boundary condition for T B , then it is possible to analyze 
an attentional region taking into account the effect of the 
non-attentional region. Specifically, the attentional 

5 region is analyzed using, as the boundary condition, a 
potential or current density on r B , which has been obtained 
by analyzing the non-attentional region. When the non- 
attentional region is further analyzed based on the 
potential or current density on the obtained division plane, 

10 the analysis of all the regions is completed. 

The object to be analyzed comprises a vertical- 
shaft pump having a diameter of 200 mm and a length of 6000 
mm, as shown in FIG. 3. As shown in FIGS. 5A and 5B, the 
pump is divided into inner pump portions 15, 16 and an outer 

15 pump portion 17. The inner surface of the pump is divided 
into a guide casing 15 of complex three-dimensional shape 
comprising a combination of complex members and having a 
helical flow path, and a column pipe 16 that can be modeled 
axially symmetrically. The three regions thus divided 

20 include a pump outer surface, a guide casing inner surface, 
and a column pipe inner surface which are handled as an 
axially symmetric open region, a three-dimensional closed 
region, and an axially symmetric closed region, 
respectively. The guide casing inner surface is divided 

25 into seven flow paths by seven helical guide vanes. Since 
these flow paths are symmetrical, one of them is taken and 
divided into three-dimensional elements. 

In order to determine the relationship between a 
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potential Vs and a current density q on boundary surfaces 
r a , r b between the pump outer surface and the column pipe 
inner surface, and the guide casing, the pump outer surface 
and the column pipe inner surface are handled as non- 
attentional regions, and subjected to a boundary element 
analysis. The former is subjected to an axially symmetric 
open-region analysis, and the latter to an axially symmetric 
closed-region analysis . 

Using the obtained relationship between the 
potential and the current density as a boundary condition, 
the guide casing (attentional region) is subjected to a 
three-dimensional closed-region analysis. Using current 
densities on the boundary surfaces T a , T b obtained from this 
analysis as boundary conditions, the pump outer surface and 
the column pipe inner surface are analyzed again, whereupon 
all analyses are completed. FIGS. 6A through 6C show the 
results of the analyses as representing potential 
distributions on the pump inner and outer surfaces. 

FIG. 6A shows the shape of the pump (segment 
positions) to be analyzed. In FIG. 6A, the horizontal axis 
represents the radial position, and the vertical axis the 
axial position. FIG. 6B shows a potential distribution on 
the pump inner surface , and FIG . 6C shows a potential 
distribution on the pump outer surface. At the positions of 
sacrificial anodes 11a, lib, 11c shown in FIG. 3, the 
potential on the pump inner surface is highly negative, and 
the potential is of about - 0.4 [V] in portions other than 
the sacrificial anodes. It can thus be seen that while the 



potential would usually be of about 0 [V] if no sacrificial 
anodes were used on a stainless pump in seawater, the 
sacrificial anodes are highly effective to prevent 
corrosion. This holds true for the pump outer surface. A 
region between T a , T b is an attentional region (guide 
casing), and other regions are non-attentional regions. On 
these division planes T a , T b , the potential distribution is 
continuous. It can be seen from the continuous potential 
distribution that continuous analytic results between 
different regions can be obtained by the related analysis 
between the different regions . 

The above embodiment represents an example of the 
present invention, and can be modified in various ways 
without departing from the scope of the invention. 

Heretofore, complex fields where two or more of 
different six regions that are modeled as two-dimensional 
objects (open and closed regions) , three-dimensional objects 
(open and closed regions), and axially symmetric objects 
(open and closed regions) are present continuously have to 
be analyzed separately. In the method according to the 
present invention, it is possible to analyze those fields in 
a related fashion. 

For example, heretofore, the inner and outer 
surfaces of a vertical-shaft pump have been analyzed 
separately though they electrochemically affect each other, 
and cannot accurately be analyzed for corrosion and 
corrosion prevention. In the method according to the 
present invention, however, it is possible to analyze the 



pump inner and outer surfaces in a related manner in a 
relatively short time. While the entire pump inner surface 
has heretofore been analyzed using a three-dimensional 
closed-region model, the method according to the present 
invention makes it possible to analyze the column pipe inner 
surface of simple shape using an axially symmetric model, 
allowing the surface to be divided into elements with ease. 

If a three-dimensional region and an axially 
symmetric region are present continuously, then it has not 
been possible to accurately determine a region that can be 
modeled axially symmetrically. According to the present 
invention, however, the regions can clearly be distinguished 
from each other as described above, and a potential 
distribution and a current density distribution can easily 
be analyzed at the boundary region, so that effective 
corrosion prevention measures can be taken. 

Industrial Applicability 

The present invention relates to a method of 
analyzing corrosion and corrosion prevention with a 
computer, and can be used to predict corrosion and corrosion 
prevention for various devices such as pumps installed in 
seawater, water, and soil. The method can be used for a 
computer simulation of a system, such as a plating facility, 
a battery cell, or an electrolysis tank, where an anode and 
a cathode are macroscopically present with an electrolyte 
interposed therebetween to form a potential field, other 
than the problems of metal corrosion and corrosion 



prevention . 
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CLAIMS 

1. A method of analyzing corrosion and corrosion 
prevention in a system where two or more same or different 
regions of six regions that are modeled as two-dimensional 
objects (open and closed regions) , three-dimensional objects 
(open and closed regions) , and axially symmetric objects 
(open and closed regions) are present continuously, 
characterized by: 

dividing the system into regions and dividing each 
of the regions into elements depending on respective models 
(two-dimensional , three-dimensional, and axially symmetric) 
thereof, and using one of the regions as an attentional 
region and another as a non-attentional region; 

giving varying values of a current density or a 
potential to the elements where current densities and 
potentials on the elements which are positionally the same 
as each other on a division plane disposed between and 
shared by two regions are equivalent to each other, the 
relationship between the current densities and potentials on 
the elements being unknown; 

solving a discrete boundary integral equation 
corresponding to the non-attentional region for the given 
values of the current density or the potential thereby to 
express the unknown relationship between the current 
densities and potentials on the division plane with the 
known relationship between current densities and potentials 
on elements in the non-attentional region other than the 



division plane; 

performing a boundary element analysis on the 
attentional region using the determined relationship between 
the current densities and potentials on the division plane 
as a boundary condition on the division plane thereby to 
determine potential and current density distributions in the 
attentional region in its entirety; and 

performing a boundary element analysis on the non- 
attentional region again using the determined potentials or 
current densities on the division plane thereby to analyze 
the regions in a related manner. 

2. A method of analyzing corrosion and corrosion 
prevention in a system, where two or more non-attentional 
regions are present continuously with respect to one said 
attentional region, characterized by: 

determining the relationship between current 
densities and potentials with respect to both division 
planes in the method according to claim 1 , 

analyzing the attentional region using the 
determined relationship as a boundary condition to determine 
potential and current density distributions in the 
attentional region in its entirety, and 

performing a boundary element analysis on the non- 
attentional region again using the determined potentials or 
current densities on the division planes thereby to analyze 
the regions in a related manner. 



3. A method of analyzing corrosion and corrosion 
prevention of a pump according to claim 2, characterized in 
that said attentional region comprises a closed region of 
three-dimensional shape in a guide casing of a pump, and the 
non-attentional regions comprise an axially symmetric closed 
region of a column pipe on an inner surface of the pump 
which is contiguous to said guide casing, and an axially 
symmetric open region of an outer surface of the pump. 

4. A method of analyzing corrosion and corrosion 
prevention of a pump according to claim 3, characterized in 
that anodes are disposed in some or all of the regions of 
the pumps and a corrosion prevention effect of the anodes is 
evaluated. 



ABSTRACT 



The present invention provides a method of 
analyzing corrosion and corrosion prevention in a system 
5 where two or more same or different regions of regions that 
are modeled as two-dimensional , three-dimensional , and 
axially symmetric regions are present continuously. 

The entire system is divided into regions, and each 
of the regions is divided into elements depending on 

10 respective models thereof, with one of the regions used as 
an attentional region and another as a non-attentional 
region. A division plane is shared by the two regions, and 
current densities and potentials on the elements which are 
positionally the same as each other on the division plane 

15 are equivalent to each other. A discrete boundary integral 
equation corresponding to the non-attentional region is 
modified to express the unknown relationship between current 
densities and potentials on the division plane with the 
known relationship between current densities and potentials 

20 on elements in the non-attentional region other than the 
division plane. A boundary element analysis is performed on 
the attentional region using the determined relationship 
between the current densities and potentials on the division 
plane as a boundary condition thereby to determine potential 

25 and current density distributions in the attentional region 
in its entirety. Then, a boundary element analysis is 
performed on the non-attentional region again using the 
determined potentials or current densities on the division 
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plane thereby to analyze the regions in a related manner. 
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Declaration and Power of Attorney For Patent Application 
Japanese Language Declaration 



^ KiTlESftaLfco 



As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated 
next to my name. 

I believe I am the original, first and sole inventor (if only one 
name is listed below) or an original, first and joint inventor (if 
plural names are listed below) of the subject matter which is 
claimed and for which a patent is sought on the invention 
entitled. 

METHOD OF ANaT.VZTNa TOPPOSTON 

AND CORROSION PREVENTION 



the specification of which 
□ is attached hereto. 

Qj was filed on S^pt^mb^r- 20. 1999 

as United States Application Number or 
PCT International Application Number 
PCT/JP99/05110 and was am ended on 

April 4, 2000 (if applicable). 

I hereby state that I have reviewed and understand the 
contents of the above identified specification, including the 
claims, as amended by any amendment referred to above. 

I acknowledge the duty to disclose information which is material 
to patentability as defined in Title 37, Code of Federal 
Regulations, Section 1 .56. 



Page 1 of 4 



8/97 



Japanese Language Declaration 



IWt, #B«**35II119* (a) - (d) BXtt365* (b) fitc 



Prior Foreign Application® 

j*B-c©$fefrUJ« 

(Number) 



(Number) 



Japan 



(Country) 
<B«) 



(Country) 
(B«) 



&f±, »35«*HSftll9* (e) S^S^TTK^H^ffr 



(Application No.) 

(ffifflS^-) 



(Filing Date) 



©65* (c) *cafr^<*?u*ccciab*ro ^ffi^s 



(Application No.) 



(Filing Date) 



(Application No.) 
(ttt*##) 



(Filing Date) 
(BJ*H) 



1001 



I hereby claim foreign priority under Title 35, United States Code, 
Section 119 (a)-(d) or 365(b) of any foreign application(s) for patent 
or inventor's certificate, or Section 365(a) of any PCT International 
application which designated at least one country other than the 
United States, listed below and have also identified below, by 
checking the box, any foreign application for patent or inventor's 
certificate, or PCT International application having a filing date 
before that of the application on which priority is claimed. 



30 /September/ 1998 
(Day/Month/Year Filed) 

(Day/Month/Year Filed) 



Priority Claimed 

a ° 

Yes No 

n □ 

Yes No 



I hereby claim the benefit under Title 35, United States Code, 
Section 119(e) of any United States provisional application(s) listed 
below. 



(Application No.) 



(Filing Date) 
(tBRH) 



I hereby claim the benefit under Title 35, United States Code, Section 
120 of any United States application®, or Section 365(c) of any PCT 
international application designating the United States, listed below 
and, insofar as the subject matter of each of the claims of this 
application is not disclosed in the prior United States or PCT 
International application in the manner provided by the first paragraph 
of Title 35, United States Code Section 112, 1 acknowledge the duty 
to disclose information which is material to patentability as defined in 
Title 37, Code of Federal Regulations, Section 1.56 which became 
available between the filing date of the prior application and the 
national or PCT International filing date of application. 



(Status: Patented, Pending, Abandoned) 



(Status: Patented, Pending, Abandoned) 

I hereby declare that ail statements made herein of my own 
knowledge are true and that all statements made on information 
and belief are believed to be true; and further that these statements 
were made with the knowledge that willful false statements and the 
like so made are punishable by fine or imprisonment, or both, under 
Section 1001 of Title 18 of the United States Code and that such 
willful false statements may jeopardize the validity of the application 
or any patent issued thereon. 
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Japanese Language Declaration 



; (B$l 



lit) 

POWER OF ATTORNEY: As a named inventor, I hereby appoint 
the following attomey(s) and/or agent(s) to prosecute this 
application and transact all business in the Patent and Trademark 
Office connected therewith: (list name and registration number) 



Norman E Obion, Reg. No. 24,61^8^ Marvin J. Spivak, Reg. No. 24,913; ; C Irvin McCieifand, Reg. No. 2 1,t24;j 3reqory J. Maier, Reg. No. 
25,599; A rthur I. Neustadt, Reg. No. 24,854jjRichard D. Kelly, Reg. No. 27,757; James D. Hamilton, Reg. No. 28,421 ; Eckhard H. Kuesters, 
Reg. No. 28,870; Robert T Pous, Reg. No ._29,099 ; Charles L. Gholz, R egHTo . J?£L395^W illiam E. Beaumont,"Reg"T5o. 30,996; Jean-Paul 
Laval leye,~Reg. No. 31 ,451; S tephen G. Baxter, Reg. No. 32,884: Richard L. Treanor, Reg. No. 36,379; Steven P. Weihrouch, Reg. No. 
32,829jj John T. Goolkasian, Reg. No. 26,142; R ichard L. Chinn, Reg. No. 3 4J305; S teven E. Lipman, Reg. No. 30,01 1 ; Carl E. Schlier, Reg. 
'4\ Np. 34.426: J ames J. Kulbaski, Reg. No . 34,648; Richard A. Neifeld, Reg. No. 35^99; J. Derek Mason, Reg. No. 35,270; Surinder Sachar, 
Iteg. No. 3 4,423; Jeffrey B. Mclntyre, Reg. No. 36,867; William T. Enos, Reg. No. 33J28; Michael E. McCabe, Jr.r^e§TN6. 37,182; Bradley 
D. Lytle, Reg. No. 40,073; and Michael R. CaseyTHegTNo. 40,294, with full powers of substitution and revocation. — — 



Send Correspondence to: 



OBLON, SPIVAK, McCLELLAND, MAIER & NEUSTADT, P.C. 



175 5 JEFFERSON DAVIS HIGHW AY 



Direct Telephone Calls to: (name and telephone number) 
(703) 413-3000 



) ~zrd 



Full name of sole or first joint inventor 

Matsuho MIYASAKA . 



Inventor's signature >7 



1 D/~in^r5V«iniVi f~*ir-^.' 



-Date 

March 26, 



2C01 



ffiJJf 



Residence c/o Ebara Research Col, Ltd., 2-1, 
Honfu j isawa 4-chome, Fujisawa-shi, JSanagasaza^ 



Citizenship 

Japan 



251-8502 Japan 



55 .*ffi©SB5fe 



Post Office Address 

same as above 



Full name of second joint inventor, if any 

Hirokazu TAKAYAMA 



Second joint Inventor's signature. 



Date 

March 26, 



20C 



mm 



Residence c/o Ebara Corporation, 11-1, 
Haneda Asahi-cho, Ohta-ku, Tokyo 144-8510 



Japan 



Citizenship 

Japan 



Post Office Address 

same as above 



(^HJW^©£tHf§$i#k:ooT fef^CcfElSLs W^ST-S Z t) (Supply similar information and signature for third and subsequent 

joint inventors.) 
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Full name of third joint inventor, if any 

Kenii AMAYA 


3: 






Third joint Inventor's signature^ A h l/\ Date 

fyffr 1 f & March 26. 20 




Residence #202 Mezon Kitto, 786 Hachiinan-cho 
Kamiinaruko, Nakahara-ku . Kawasaki -shi r K^n^rj 


mm 


Citizenship 211-0001 Japan _ 0 
Japan ""J 








Post Office Address 
same as above 









Full name of fourth joint inventor, if any 

Shiaeru AOKI 




Fourth joint Inventor's signature jiu ^ Date 

ft^-^X- March 26, 20 




Residence 2-4-302, Wakabadai, Asahi-ku, 
Yokohama-shi, Kanaqawa 241-0801 Japan 




Citizenship 
Japan 




Post Office Address 
same as above 










Full name of fifth joint inventor, if any 




Fifth joint Inventor's signature Date 




Residence 




Citizenship 




Post Office Address 










Full name of sixth joint inventor, if any 




Sixth joint Inventor's signature Date 




Residence 


an 


Citizenship 




Post Office Address 






(SB A i tz ti-f tifeittO*|pl5BW#t-J* L.T (Supply similar information and signature for third and sub- 
SB*? <£ € if f#~f ZZto) sequent joint inventors.) 
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